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Abstract 

Magnetosonic waves propagating perpendicular to an external magnetic field are studied 
with attention to the effect of multiple ion species. First, power spectra of magnetic field fluctu- 
ations and autocorrelation functions in thermal equilibrium plasmas are numerically obtained. 
In a multi-ion-species plasma, besides u ~ kvA mode, numerous waves are present near many 
different ion cyclotron frequencies. The autocorrelation function of the quasi-mode consisting of 
these waves is not recovered to its initial value, owing to the phase mixing of these waves. Next, 
with particle simulations, evolution of a macroscopic perpendicular disturbance is investigated. 
In a multi-ion-species plasma, this disturbance is damped. The energy is transferred to from 
the magnetic field to the ions. 
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1. Introduction 

The presence of multiple ion species introduces many interesting effects on magnetosonic 
waves [1-7]. For instance, in a two-ion-species plasma, the magnetosonic wave is split into two 
modes. Nonlinear pulses of these modes are damped, even when they propagate perpendicular 
to the magnetic field [5-7]. The damping is due to energy transfer from the pulse to heavy ions 
[8,9]. Periodic waves are not damped even in this case. However, the collective behavior of these 
waves in a multi-ion-species plasma would be different from that in a single-ion-species plasma. 

Recently, a study has been made on collective behavior of ion Bernstein waves in thermal- 
equilibrium plasmas with multiple ion species [10]. Each perpendicular ion Bernstein wave with 
u ~ nil; is undamped in a collisionless plasma [11], where f2j is the ion cyclotron frequency and 
n is the integer. In a single- ion-species plasma, the autocorrelation function of the quasi-mode 
consisting of these waves shows periodic behavior with time period 27r/J7j. On the other hand, in 
a multi-ion-species plasma, the autocorrelation function is initially damped and is not recovered. 
This is caused by the phase mixing of numerous waves excited at the harmonics of many different 
ion cyclotron frequencies. This damping mechanism could be important in space plasmas where 
many ion species exist, with each species having many different ionic charge states. 

In this paper, we study perpendicular magnetosonic waves in thermal-equilibrium, multi-ion- 
species plasmas where each particle species has its own Maxwellian velocity distribution. We 
assume that all the ion species have an equal temperature, while electrons can have a different 
temperature because relaxation time between electrons and ions via collisions is very long. 
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In Sec. 2, we numerically calculate power spectra and autocorrelation functions of magnetic 
field fluctuations due to the magnetosonic waves. In a single-ion-species plasma, the autocorre- 
lation function is not damped, because the wave with uo ~ kv\ is dominant mode. Here, va is 
the Alfven speed and k is the perpendicular wavenumber. On the other hand, in a multi-ion- 
species plasma, besides this mode, numerous waves are present near many different ion cyclotron 
frequencies. Owing to the phase mixing of these waves, the autocorrelation function does not 
return its initial value. In Sec. 3, evolution of a macroscopic disturbance and associated energy 
transport are studied by particle simulations. In a multi-ion-species plasma, the macroscopic 
disturbance is damped, and the energy is transferred from the magnetic field to the ions. 



2. Numerical Calculation 

We consider extraordinary waves propagating perpendicular to an external magnetic field with 
frequensies smaller than lower hybrid frequency; we call these waves as magnetosonic waves. The 
dispersion relations of the magnetosonic waves are given by 

Ans = e 2 xy /e xx + e yy - c 2 k 2 /uo 2 = 0, (1) 

where c is the light speed, and e xx ,e xy , and e yy are defined as 
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Here, the subscript j refers to electrons (e) or ion species (H, He, C, • • •), £lj is the cyclotron 
frequency, and uo p j is the plasma frequency. Also, T n (pj) = I n (fij) exp(—/j,j), where /„ is the 
modified Bessel function of the nth order, and fij = k 2 with pj the gyro-radius. 

The fluctuation spectrum of magnetic fields due to the magnetosonic waves in a spatially 
homogeneous, thermal equilibrium plasma is written as 
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where fee is the Boltzmann constant, and uj n is the roots of the dispersion relation D ms =0. 

In a multi-ion-species plasma with H being major ions, there are three kinds of waves in the 
long wavelength region, pi <C 1; 

uj ~ kv\, (7) 



oj 2 M-u 

to ~ n s + (n H - n s ), (8) 

and 

uo ~ nfli. (9) 
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Here, the subscript s in Eq. (JHJ) refers to heavy ion species (He, C, O, • • •). The waves with Eqs. 
l [7| ) and (JHJ) exist even in a cold plasma [4] , while the mode with Eq. (JHJ) is caused by ion kinetic 
effects [12]. 

For a given wavenumber k, there are many waves with different frequencies. Autocorrelation 
function of the quasi-mode consisting of these waves is obtained from P[ui) though the Fourier 
transformation in co as 

/■oo 

Ck{r) = \ \B k ,uj\ 2 exp(-iuT)du. 



(10) 



We pay attention to how values of |Cfc(r)| are reduced by the presence of multiple ion species. 
The reduction of |C&(t)| indicates that energy transport can be enhanced [13]. 

We numerically calculate specific values of P(u>) and Cfc(r) for three different plasmas and 
compare them. [In the calculation, we retain the terms from n = —10 to 10 for the ions and 
the n = and 1 terms for the electrons in Eqs. (J2J), ([3*]). and ©.] The three plasmas that we 
consider are single-ion (H + ), three ion (H + , He +2 , C +5 ) and six-ion (H + , He +2 , C +6 , +6 , Si +9 , 
and Fe+ 13 ) species plasmas. The cyclotron frequencies of these ions normalized to Oh are taken 
to be He = 0.5, n c = 0.417, o = 0.375, Si = 0.321, and Fe = 0.232. The densities of the 
ions normalized to ran are nHe = 0.1, nc = no = 0.01, and rasi = np e = 0.005. The magnetic 
field strength is |f2 e |/w P e = 1. The plasma beta value is /3 = 0.0625. 
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Figure 1: Power spectra of magnetic field fluctuations with kpn = 0.1 in three different 
plasmas. 

Figure 1 shows power spectra of the mode with kpn = 0.1 in the single-, three, and six- 
ion-species plasmas, where P(uj) is normalized to 7t/cbT '. In the single-ion-species plasma, the 
wave with u ~ kv^ (~ 0.42Oh) is the dominant mode. Even though there are waves with 
ijj ~ nf^Hj their amplitudes are quite small. In the three- ion-species plasma, besides the u ~ kvA 
(~ 0.39^h) mode, the waves near ^He and fic are present; these frequencies are given by Eq. (JHJ 
with s = He or C. The amplitudes of the waves with u ~ n^H) n ^Re, and n^c are much smaller. 
In the six-ion-species plasma, the waves near fio, ^Si> an d also exist. Their amplitudes are 
not small, although the abundances of the heavy ions are very small. 
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Figure 2: Autocorrelation functions of the fluctuations with kpn = 0.1 in the same plasmas 
as shown in Fig. 1 

Figure 2 shows autocorrelation functions normalized to their initial values Cfc(0). In the 
single-ion-species plasma, Ck(r) oscillates with the period 2Tr/(kv\) and is undamped. In the 
three- and six-ion-species plasmas, Ck(r)'s do not return to their initial values till the end of the 
calculation. As the number of ion species increases, the amplitude of the oscillation decreases 
more quickly, owing to the phase mixing of more waves. 

If effects of collisions are entirely neglected, Ck(r) returns to its initial value at the time of 
the least common multiple of all the wave periods. However, this time is extremely long in space 
plasmas where the number of ion species is very large (moreover, each ion species has many 
different ionic charge states) and numerous waves exist. On such a long time scale, the effects of 
collisions must be important, which reduces |Cfc(r)|. Accordingly, Cfc(r) would not be recovered 
and would keep smaller values than in a single-ion-species plasma. 

3. Particle Simulations 

By means of a one-dimensional (one space and three velocity components), electromagnetic 
particle code with full ion and electron dynamics, we study collective behavior of magnetosonic 
waves in a multi-ion-species plasma. The system size is L x = 512A 9 , where A g is the grid 
spacing and is equal to the electron Debye length. We use periodic boundary conditions. The 
external magnetic field is in the z direction, and its strength is |O e |/<-^pe = 4.0. The total number 
of electrons is N c = 262, 144. The plasma (3 value is (3 = 0.03. 

We simulate single-ion (a) and four-ion (a, b, c, and d) species plasmas. We choose the mass 
ratios as m a /m c = 50, m^/m^ = m c /m a = 2, and md/m a = \/5. In order to see the effect 
of multiple ion species with a small number of ion species, we have taken the irrational ion mass 
rations for b and d ions. The charges are the same, q a = qb = Qc = Qd = \Qc\- The ion densities 
are set to be rib = n c = n d = 0.2n a . 

Firstly, we observed that autocorrelation functions of fluctuations propagatig perpendicular 
to the external magnetic field are not recovered in the four-ion-species plasma. Next, as a initial 
condition, we set the magnetic field to have a finite amplitude disturbance with a monochro- 
matic cosine profile, 5B z (x)/ Bq = 0.02 cos(fcox), where Bq is the external magnetic field and 
koPa = 0.01(fco^A = 0.68$7 a ). We then study evolution of its disturbance and associated energy 
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transport. Initially, all the ion species have equal temperature; the electron-to-ion temperature 
ratio is chosen to be Tj/T e = 0.1. 




Figure 3: Time variations of total magnetic-field energies and ion kinetic energies. The 
thin and thick lines represent single- and four-ion-species plasmas, respectively. 

Figure 3 shows time variations of total magnetic- field energy Eb and ion energy K — Kq, 
where K is the total energy of all the ions and Kq is the initial one. The energies are normalized 
to m c v\ c . The thin and thick lines denote energies in the single- and four-ion-species plasmas, 
respectively. In the single-ion-species plasma, the magnetic field energy and ion kinetic energy 
oscillate with period tt/^qva)- On the other hand, in the four-ion-species plasma, the magnetic 
field energy is reduced and does not return to its initial value. The ion energy is rapidly increased 
and then keeps large values. 

Figure 4 shows long time variations of total magnetic field energies in the single- and four-ion- 
species plasmas. (For the single-ion-species plasma, the values averaged over the time period 
t = 27r/(fcofA) are plotted.) The magnetic field energy in the single-ion-species plasma keeps 
almost constant even for the long period. On the other hand, in the four-ion-species plasma, 
the magnetic-field energy is eventually damped. 

Figure 5 shows time variations of ion kinetic energies, which is defined as 

KE = J2 J dx J dvmifi(x,v)(v- < v^x) >) 2 , (11) 

i 

with < Vi(x) > the fluid velocity at position x, and KE(0) is the initial value of KE. In 
the four-ion-species plasma, the ion kinetic energy is increased. Evidently, the presence of the 
multiple ion species can enhance the energy dissipation. 

4. Summary 

We have studied collective behavior of perpendicular magnetosonic waves in multi-ion-species 
plasmas. We have numerically shown that the autocorrelation functions in a thermal-equilibrium 
plasma are not recovered, because, in addition to the to ~ kvA mode, many waves exist near 
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Figure 4: Long time variations of total magnetic-field energies. The solid line and dashed 
lines represent the energies in the four- and single-ion-species plasmas, respectively. For 
the single- ion-species plasma, the values averaged over the time period t = 2n/k VA are 
plotted. 
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Figure 5: Long time variations of ion kinetic energies in the single- and four- ion-species 
plasmas. 
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many different cyclotron frequencies. Furthermore, we have shown with particle simulations, 
that the macroscopic disturbance is also damped in a multi-ion-species plasma, and that the 
energy is transferred from the magnetic field to the ions. We have not yet understood the transfer 
mechanism. To do this, we will further investigate with particle simulations, for example, how 
the energies depend on the initial conditions or on the ion species. 
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